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ABSTRACT 

In  recent  years,  a  compressive  sensing  based  underwater  imaging  system  has  been  under  investigation:  the  Compressive 
Line  Sensing  (CLS)  imaging  system.  In  the  CLS  system,  each  line  segment  is  sensed  independently;  with  regard  to 
signal  reconstruction,  the  correlation  among  the  adjacent  lines  is  exploited  via  the  joint  sparsity  in  the  distributed 
compressive  sensing  model.  Interestingly,  the  dynamic  compressive  sensing  signal  model  is  also  capable  of  exploiting 
the  correlated  nature  of  the  adjacent  lines  through  a  Bayesian  framework.  This  paper  proposes  a  new  CLS  reconstruction 
technique  through  the  integration  of  these  different  models,  and  includes  an  evaluation  of  the  proposed  technique  using 
the  experiment  dataset  obtained  from  an  underwater  imaging  test  setup. 


1.  BACKGROUND 

The  laser  line  scanner  (LLS)  is  state-of-art  underwater  active  imaging  technology  [6].  The  LLS  system  employs  serial 
raster  scan  image  formation  and  is  compatible  with  the  push-broom  type  image  acquisition  commonly  used  in  many 
survey  operations.  While  the  LLS  system  provides  high  photon  efficiency  and  narrow  instantaneous  field  of  view  that 
are  critical  for  operating  in  degraded  visual  environment,  the  non-adaptive  serial  raster  scan  image  formation  can  be  a 
concern  in  certain  situations.  For  example,  higher  repetition  laser  will  be  required  with  increased  platform  speed  to 
maintain  the  same  image  spatial  resolution,  which  inevitably  will  increase  the  system  cost  and  require  more  sophisticated 
noise  mitigation.  For  unmanned  platforms  such  as  UAVs  and  AUVs  where  power  is  at  a  premium,  another  concern  is 
that  the  LLS  sensor  may  consume  significant  power  acquiring  redundant  data  that  will  be  discarded  during  the  image 
compression  phase  (Figure  1).  _ 


Reconstructed 


— *  Output 


Figure  1.  Traditional  image  acquisition  paradigm  in  a  LLS  system  [3] 


The  compressive  line  sensing  (CLS)  imaging  system  has  been  developed  [1]  to  address  these  concerns.  By  adopting  the 
concept  of  resource  compression,  this  system  aims  to  achieve  faster  signal  formation,  more  flexible  system  configuration 
and  higher  energy  efficiency.  The  technique  also  enables  compact  and  robust  system  design.  A  series  of  test-tank 
experiments  were  conducted  to  validate  the  CLS  concept  [3],  and  while  the  results  demonstrated  the  effectiveness  of  the 
underlying  distributed  compressive  sensing  (DCS)  signal  model  was  demonstrated  in  [3],  there  was  still  room  for 
improvements  in  signal  reconstruction  process.  In  this  paper  we  describe  our  attempt  to  improve  the  performance  of  the 
CLS  system  by  integrating  the  “dynamic”  compressive  sensing  concept  into  the  signal  reconstruction. 


We  provide  some  theoretical  foundations  in  Section  2  and  then  introduce  the  CLS  concept  in  Section  3.  In  Section  4  we 
introduce  the  new  signal  reconstruction  technique  through  the  integration  of  distributed  compressive  sensing  and 
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dynamic  compressive  sensing  signal  models.  In  Section  5  we  apply  this  technique  to  some  test-tank  datasets  and  provide 
results  and  detailed  analysis.  Our  conclusions  are  summarized  in  Section  6. 


2.  CS  AND  DCS  THEORIES 

Compressive  sensing  (CS)  is  a  framework  for  the  simultaneous  sampling  and  compression  of  sparse  signals  using 

incomplete  linear  measurements  [5,  9].  A  iC-sparse  signal  X  G  RNxl  (i.e.,  the  projection  of  X  on  a  sparsifying  basis  T7  G 
Rnxn:X  =  x¥a  contains  K 

«  N  non-zero  entries)  can  be  recovered  with  overwhelming  probability  using  more  than 
M  =  0(K\ogN)  incoherent  linear  measurements:  y  =  AX  =  AT7# ,  when  the  measurement  matrix  A  G  RMxN  is 
incoherent  with  the  sparsifying  basis  T7  and  satisfies  the  Restricted  Isometry  Property  (RIP)  [5].  The  vector  a  (therefore 
X)  can  be  recovered  from  the  measurements  y  by  solving  an  L 1 -minimization  problem.  Random  bases  such  as  Gaussian 
and  Bernoulli  are  common  measurement  matrices  that  satisfy  the  RIP  constraint.  For  natural  signals,  the  Fourier,  discrete 
cosine  transform  or  wavelet  domains  are  sparsifying  bases  commonly  used  in  various  image/video  compression 
standards  and  codecs. 

The  CS  theory  essentially  exploits  the  intra-signal  redundancy  within  a  single  source,  and,  there  has  been  significant 
interest  in  extending  it  to  cope  with  the  inter-source  correlations.  We  briefly  introduce  two  such  approaches  herein.  One 
idea  is  based  on  the  distributed  compressive  sensing  concept,  which  we  abbreviate  as  DCSa.  The  other  school  of  thoughts 
includes  modified  compressive  sensing  residual  (Lu  et  al.  [14])  and  dynamic  compressive  sensing  (Ziniel  and  Schniter 
[18]),  we  abbreviate  these  techniques  as  DCSb. 

DCSa  is  closely  related  to  the  distributed  source  coding  theorems  (Slepian  and  Wolf  [15]  and  Wyner  and  Ziv  [17]), 
which  hold  that  the  minimum  rate  of  encoding  statistically  dependent  sources  independently  is  the  same  as  the  minimum 
rate  of  encoding  them  jointly,  when  the  sources  are  decoded  jointly  and  their  differences  are  Gaussian. 


DCSa  attempts  to  exploit  the  inter-signal  redundancy  among  distributed  and  correlated  sources  through  the  establishment 
of  the  proper  joint  sparsity  models  (JSMs)  [4].  In  JSM-1  is  framework  adopted  in  the  CLS  framework,  all  sources 
Xt  within  the  group  consist  of  a  sum  of  a  common  component  Zc  and  a  component  that  is  unique  to  each  source  Zl:Xl  = 
Zc  +  Zlf  l  =  0,1 ...  L,  where  L  is  the  number  of  signals  to  be  solved  jointly.  On  a  sparsifying  basis  T7,  Zc  and  Zt  can  be 
expressed  as  Zc  =  T7^,  ||crc||0  =  Kc ,  and  Zt  =  T7^,  ||aq||0  =  Kt  and  both  are  sparse,  i.e.,  Kc  «  N  and  Kt  «  N.  The 
matrix  form  of  the  measurement  model  is  illustrated  below: 
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where  y  =  [y(  1), ,  y(L)]r  G  RLMxlare  the  measurements  of  L  sources  within  the  group;  a  =  [a(l), ... ,  a(L )]r  G 
rLNx1  are  the  coefficients  of  the  sources  on  the  sparsifying  basis  T7,  and  A  is  the  measurement  matrix.  Solving  such 
JSM-1  problem  can  be  casted  as  LI  minimization  with  the  cost  function  formulated  by  revising  Equation  (1)  [4]: 


a*  =  arg min  ||ac| |x  +  ||a1||1  +  ---+  HaJi 
subject  to  ||y  -  A^al^  <  £ 


where  \\ai\\1  =Z?Lil&z(0l  is  the  LI -norm  of  aL. 

DCSb  formulates  a  group  of  correlated  sources  as  recursively  reconstructing  time  sequences  of  sparse  spatial  signals. 
With  regard  to  the  measurement  acquisition,  DCSa  and  DCSb  are  identical  -  in  that  each  source  is  measured 
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independently.  With  regard  to  the  signal  reconstruction,  both  DCSa  and  DCSb  assume  there  is  a  common  component 
among  the  different  sources.  The  key  difference  is  that  DCSb  uses  the  knowledge  of  the  source  at  the  previous  time 
instant  as  a-priori  to  increase  the  sparsity  of  the  signal  at  current  time  instant.  As  a  result,  the  signal  reconstruction 
performance  is  improved. 

There  are  two  approaches  in  DCSb  -  modified-CS-residual  [14]  and  dynamic  CS  [18].  In  modified-CS-residual,  the 
support  estimate  { S }  from  the  previous  time  instant  will  be  used  as  the  known  part  of  the  support.  Then  at  the  current 
time  instant,  the  scheme  attempts  to  solve  the  signal  that  satisfies  the  data  constraint  and  whose  support  contains  the 
smallest  number  of  new  additions  to  the  known  support  {' T }  [14].  The  scheme  can  be  described  by  Equation  (3)  [14]: 

a;  =  argmin  rUCa^c  J^  +  ||ytres  -  A^a||2 

—  & t,temp  dtp 

where  ytres  =yt-  &Wat>temp  ,  attemp  is  constructed  using  support  from  step  t-1  St_p.  at  temp  =  (at_1)5t_1. 

In  Ziniel  and  Schniter  [18],  the  dynamic  CS  adopts  a  slight  different  approach.  One  of  the  main  motivations  was  to  tackle 
the  Bayesian  inference  high-dimensional  problems  in  a  computationally  efficient  manner.  A  probabilistic  dynamic  CS 
signal  model  that  captures  both  amplitude  and  support  correlation  structure  was  proposed: 

a?  =  (4) 

where  af  represents  nth  coefficient  at  time  instant  t,  {s}^=1  and  {0}f=1  are  two  hidden  random  processes  and  T  is  the 
duration  of  the  time  sequence.  The  binary  vector  {s}  describes  the  support  of  {a}.  {0}  follows  normal  distribution  to 
describe  the  amplitudes  of  the  active  elements  of  {a}.  One  critical  aspect  of  this  algorithm  is  to  update  { at  }  and  to 
propagate  the  updates  to  the  next  time  instant  efficiently.  A  novel  approximate  message  passing  (AMP)  algorithm  -  an 
unconventional  form  of  loopy  belief  propagation,  was  adopted  to  perform  soft  signal  estimation  and  support  detection 
with  a  computational  complexity  that  is  linear  in  all  problem  dimensions  [18]. 


3.  CLS  IMAGING  SYSTEM  THROUGH  SCATTERING  MEDIUM 
3.1.  CLS  imaging  system  architecture 

The  CLS  imaging  system  can  be  conceptualized  as  a  “coding”  machine  (Figure  2).  During  image  acquisition,  entries 
from  a  codebook  {A}  (i.e.,  a  measurement  pattern  A*)  sequentially  encode  the  target.  During  image  reconstruction,  the 
encoder  output  {Y}  and  an  appropriate  codebook  {®}  can  used  to  recover  the  target  scene.  Here  the  purposely  different 
notation  ®  indicates  that  the  encoder  and  decoder  may  require  different  codebooks  when  imaging  through  the  scattering 
medium. 


Reconstructed 
_ I  image 

Model 

Prediction 

I  DECODER 

Decoding  — 


Figure  2.  CLS  System  as  a  Coding  Machine 

The  CLS  illuminator  shares  some  similarity  with  the  streak  tube  imaging  Lidar  [15].  Nevertheless,  instead  of  one¬ 
dimensional  (ID)  fan  beam,  a  laser  source  modulates  the  digital  micromirror  device  DMD  and  projects  a  series  of  ID 
“fan  ’’patterns  onto  the  current  target  line  (Figure  3).  Each  pattern  corresponds  to  one  row  in  the  CS  measurement  matrix 
(i.e.,  codebook).  At  the  receiver,  similar  to  the  LLS  system,  a  single  element  receiver  such  as  the  photomultiplier  tube 
(PMT)  records  the  total  photon  reflection  corresponding  to  the  modulation  of  the  spatial  pattern  with  the  target  as  the 
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measurement.  Each  line  is  measured  independently. ;  the  platform  motion  enables  the  sensing  of  the  target  scene  in  line- 
by-line  fashion. 

In  the  CLS  system,  resource  compression  is  achieved  by  reducing  the  required  measurements  per  line  while  still  maintain 
desired  image  quality. 
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3.2.  CLS  Signal  Model 


Figure  3.  CLS  Imaging  system  architecture 


For  the  CLS  system  operating  in  the  scattering  medium,  such  as  in  the  underwater  environment,  when  a  binary  pattern  of 
highly  collimated  light  is  emitted  from  the  transmitter,  the  light  will  diffuse  and  attenuate  as  it  propagates  to  the  target. 
The  degree  of  spreading/attenuation  is  determined  by  the  laser  beam  divergence,  target  range,  and  the  inherent  optical 
properties  of  the  water.  Assuming  a  wide  receiver  aperture  and  Lambertian  reflection  at  the  scene,  the  total  photon  flux 
corresponding  to  the  mth  measurement  Y  can  be  represented  by  the  equation: 


Ym(r  +  A  rm)  =  Ps 


sys 


(BSFIT(i,j)  *  Am(i,j ))  °X(i,r  +  A rm 


+ 


(5) 


Where  *  denotes  convolution  and  °  represents  the  Hadamard  point-by-point  product,  Ym  is  the  total  reflected  photon 
flux;  r  is  the  current  line  under  investigation;  Psys  is  a  constant  accounting  for  the  system  configuration  (receiver 
aperture,  illumination  power  etc.);  aTR  is  an  attenuation  coefficient  related  to  radiative  transfer  from  the  scene  back  to  the 
sensor;  =  Am  *  BSFIT  is  the  measurement  pattern  on  the  target  (i.e.,  the  original  binary  pattern  after  propagating 
from  the  illuminator  to  the  target  through  a  scattering  and  absorbing  medium  described  by  the  beam  spread  function 
BSFjt );  Om  o  X  is  the  (information  bearing)  “imprinted”  pattern  due  to  the  modulation  of  the  measurement  pattern  by  the 
target  scene  reflectance  pattern;  N  is  the  number  of  pixels  within  each  line;  N^m  is  the  number  of  ‘on’  pixels  in  the 
illumination  pattern  m;  U  is  the  effective  vertical  aperture  due  to  the  spreading  in  the  vertical  (along  track)  direction; 
A Lm  =  mvAt  is  the  platform  displacement  when  the  measurement  pattern  m  is  projected  (relative  to  the  previous 
pattern);  At  is  the  refresh  time  of  the  spatial  light  modulator  (SLM)  device  and  v  is  the  platform  speed. 

We  will  make  two  reasonable  approximations:  1)  the  SLM  refresh  rate  is  sufficient  fast  relative  to  the  platform  speed; 
and  2)  the  beam  spread  function  has  a  separable  kernel  that  can  be  decomposed  into:  BSFIT  «  BHBV ,  where  Bv  and  BH 
describe  the  vertical  and  horizontal  beam  spreading  respectively.  Taking  into  consideration  that  Am  is  now  a  ID  pattern, 
Equation  (5)  can  be  simplified  to: 
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4.  NEW  CLS  SOLVER  BASED  ON  THE  INTEGRATION  OF  DCSa  AND  DCSb 
4.1.  DCSa  JSM-1  based  CLS  Solver 


jointly.  The  corresponding  CS  cost  function  then  becomes: 

a*  =  argmin  HaJi  +  ||a1||1  +  -+  HaJi 
subject  to  ||  V  —  $lTa||2  <  £,  X  =  '¥a>  0 

where  L  is  the  group  of  lines  to  be  solved  jointly.  Compared  to  Equation  (1), 
Equation  (7)  also  incorporates  a  non-negative  constraint.  In  order  to  comply 
with  the  DCSa  measurement  matrix  formation  in  Equation  (1),  we  introduce 
another  approximation,  ignoring  the  vertical  spreading: 

Am(r)  (8) 

The  corresponding  process  is  illustrated  in  the  flow  chart  in  Figure  4. 
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Figure  4.  CLS  System  processing  flow 


One  unique  feature  of  this  framework  is  that  it  shares  some  similarity  with 
the  “lucky  imaging”  scheme.  As  illustrated  in  Equation  (9),  assuming  the 

current  line  group  parameter  L  >  1,  then  the  reconstruction  of  line  r  will  be  included  in  L  solutions: 


a\  —  argmin 


Kill  +  IK-rlll  - f 


a?  =  argmin  ||a? IU  +  ||ai,_T-+1  IU  +  • 


||ar||1(  subject  to  \\Y[_r -$^a\\2  <£,X1  =  ^a1  >  0 

|ttr||  “I-  ||&r+i Hi, 


subject  to  HT/V+i  -  $  Vallz  <  £, X2  =  'Ta2  >  0 


(9) 


at 


argmin  Htfflh  +  Hor-lli  H — l|arll  +  llar+illi<  subject  to  ||Trr+i  -0  'T'aMz  <  £,XL  =  ^VaL  >  0 


We  will  buffer  these  L  solutions  and  the  optimum  reconstruction  for  line  r  can  then  be  derived  from  these  L  solutions. 
One  straightforward  approach  is  to  apply  filtering  such  as  a  median  filter  to  obtain  the  final  solution: 

a£  =  median[ar],  foru  =  1  ...L 


While  this  framework  has  been  validated  through  application  to  test-tank  experiment  datasets  acquired  at  different 
turbidity  cycles  (Figure  5)  [3],  the  image  quality  at  higher  turbidities  (Figure  5b  and  5c)  leaves  much  room  for 
improvement. 


(a)  Clearwater  (b)  C=0.6  (3  attenuation  length)  (c)  C=0.8  (4  attenuation  length) 

Figure  5.  Reconstructed  image  at  2:1  compression  using  test  tank  experimental  data  (Clearwater  dataset  acquired  on  9/20/2013,  two 
turbidity  cycle  datasets  acquired  on  4/29/2014.  Two  different  optical  front  end  configurations  were  used  in  these  two  experiments)  [3] 


The  deficiency  in  the  illuminator  optical  front  end  of  the  prototype  system  is  one  of  two  factors  impacting  image  quality 
[3].  On  the  algorithm  side,  as  shown  in  Equation  (6),  the  measurements  of  every  line  are  affected  by  the  adjacent  lines 
and  pixels  (i.e.,  forward  scattering),  and  the  non-information  bearing  backscattering  photons.  While  the  forward 
scattering  among  the  adjacent  pixels  within  the  same  line  were  addressed  in  Equation  (7),  the  contribution  from  the 
adjacent  lines  is  problematic.  The  forward  scattering  from  the  adjacent  lines  essentially  means  that  the  underlying 
process  is  non-causal  since  the  measurements  for  the  current  line  contain  information  from  the  previous  lines  (i.e.,  past) 
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and  the  subsequent  lines  (i.e.,  future).  The  DCSa  framework  (Equation  1)  describes  a  memoryless  system,  which  is  the 
need  for  the  approximation  in  Equation  (8).  With  regard  to  backscattering  mitigation,  although  the  polarity  flipping 
technique  was  shown  to  be  helpful  [1],  the  impact  from  the  backscattering  remained  prominent,  in  Figure  5b  and  5c. 

While  there  certainly  is  a  need  to  improve  the  prototype  system  hardware,  we  will  focus  on  enhancing  image  quality  by 
improving  the  algorithm. 

4.2.  Proposed  CLS  Reconstruction  Technique  through  the  Integration  of  DCSa  and  DCSb 

From  the  compressive  sensing  perspective,  the  effect  of  the  backscattering 
and  forward  scattering  interference  is  reduced  signal  sparsity  via  more 
information  injected  the  measurements.  On  the  other  hand,  the  main  theme  in 
the  DCSb  school  of  thought  is  to  use  a  priori  information  to  improve  signal 
sparsity  when  recursively  solving  a  time  sequence  of  signals.  Additionally, 
both  causal  filtering  and  non-causal  smoothing  can  be  easily  realized  in 
DCSb,  overcoming  one  limitation  in  the  current  CLS  solver.  Therefore,  the 
proposed  new  design  attempts  to  better  mitigate  the  interferences  from 
backscattering  and  forward  scattering  through  the  integration  of  the  DCSb 
concept  into  the  current  DCSa  based  framework. 

The  flow  chart  in  Figure  6  outlines  the  proposed  two-stepped  solution.  The 
new  technique  consists  of  three  phases.  We  start  with  the  reasonable 
assumption  that  the  imaging  system  will  first  traverse  through  regions  with 
no  target  presence  to  perform  the  initialization,  allowing  us  to  measure  the 
contribution  of  the  volume  backscatter  /?  .  We  assume  the  volume 
backscatter  to  be  invariant  with  the  same  system  configuration  and 
environment  conditions  (i.e.,  same  IOPs).  We  then  proceed  to  invoke  the 
Level-one-solver: 

5*  =  argmin  HadU  +  ||a1||1  +  -+  HadU 
subject  to  \\Y  -&'¥a-Y$P\\2  <  £,  X  =  >  0  v±1' 

where  y  is  a  weighting  constant.  Comparing  Equations  (7)  and  (11),  this  step  is  essentially  the  original  DCSa  based 
solver  adjusted  for  the  contribution  of  the  backscattering  assessed  prior  to  the  measurements.  The  more  significant  aspect 
of  the  proposed  technique  is  that  it  also  incorporates  the  Level-two-solver,  which  becomes  involved  after  the  current  line 
is  outside  the  effective  vertical  aperture.  Consequently,  the  solution  incorporates  contributions  to  the  measurements  from 
the  adjacent  lines  obtained  in  Level-one  in  addition  to  the  backscatter: 

a;  =  argmin  llaji  +  ||a1||1  +  -+  ||aL||i 

subject  to  || Y  -dVa-  y$P  -  n[$'V(alr1_u  +  -  +  a*.1  +  -  +  a^+u)]  ||2  <  £,  X  =  'Va>  0  ^ 

a,  =  [lUr1  +  ap 

where  y  and  j.i  are  two  weighting  constants,  a*1,  i  =  r  —  U  ... ,r, ...  r  +  U  are  the  solutions  from  the  Level-one-solver. 
Equation  12  shares  strong  similarity  with  the  modified-CS-residual  in  Equation  3.  The  two  main  differences  are:l)  The 
solver  is  based  on  DCSa  and  2)  In  the  final  solution,  we  annihilate  the  terms  corresponding  to  the  backscatter  and 
forward  scattering  from  adjacent  lines  estimated  in  the  Level-one-solver. 

5.  EXPERIMENTAL  RESULTS 

We  now  present  the  results  of  applying  the  new  algorithm  to  the  same  dataset  that  produced  the  images  in  Figure  5.  We 
first  briefly  describe  the  prototype  systems  and  the  experimental  setup.  [3]. 


Figure  6.  Flow  chart  illustrates  the  proposed 
CLS  reconstruction  algorithm 
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5.1.  Prototype  System  Architecture  and  Experimental  Setup 

The  overall  system  consists  of  the  illumination  and  the  receiver  subsystems  (Figure  7a).  The  configuration  of  the  receiver 
consists  of  a  Hamamatsu  R9880U-210  PMT  with  12  degree  field  of  view.  The  core  of  the  illumination  subsystem  is  a 
compact  DLP  Lightcrafter  evaluation  module,  consisting  of  a  DMD  with  an  array  of  608  x  684  pixels. 

The  validation  tests  were  conducted  in  the  FAU  Harbor  Branch  optical  test  tank  (Figure  7b).  A  linear  drive  mounted 
against  the  catwalk  provides  both  vertical  and  lateral  translation.  The  linear  drive  motion  accuracy  is  +/-0.5  mm.  The 
catwalk  itself  can  be  manually  moved  to  provide  an  additional  degree  of  freedom.  During  the  experiment  a  75cm2 
USAF1952-B  target  (Figure  6d)  was  attached  to  the  linear  drive  and  placed  5m  away  from  the  view  port  (Figure  7c).  A 
Matlab  program  was  developed  to  synchronize  the  linear  drive  movement,  the  data  acquisition  and  the  pattern  cycling  on 
the  Lightcrafter.  During  the  experiments,  the  linear  drive  moved  at  1.5mm  increment.  The  NI-6133  data  acquisition 
board  operated  at  80KHz  sampling  rate.  Each  data  acquisition  was  triggered  by  the  Lightcrafter  trigger  output.  Four 
hundred  points  were  acquired  for  each  pattern,  and  the  average  of  these  400  points  was  recorded  as  the  corresponding 
measurement.  While  the  Lightcrafter  is  capable  of  a  4KHz  refresh  rate  in  the  binary  “pattern  sequencing”  mode,  the 
Lightcrafter  was  operated  at  a  400Hz  refresh  rate  to  ensure  sync  stability  during  all  the  tests. 


Figure  7  Illustration  of  the  prototype  system,  the  FAU  Harbor  Branch  Optical  Imaging  Test  Tank  and  the  test  setup  configuration 

For  image  quality  evaluation,  we  will  adopt  the  anisotropy  based  blind  image  quality  assessment  metric:  Anisotropic 
Quality  Index  (AQI)  proposed  in  [19].  AQI  computes  the  variance  of  the  expected  entropy  of  a  given  image  as  a  function 
of  the  directionality,  which  is  taken  as  an  anisotropy  indicator.  Experiments  results  in  [19]  show  this  metric  presents 
some  desirable  features  that  resemble  those  from  an  ideal  image  quality  function,  constituting  a  suitable  quality  index  to 
assess  both  the  contrast  and  resolution  of  natural  images.  AQI  produces  a  score  between  0  and  1,  with  1  being  the  higher 
score  indicates  better  image  quality/fidelity. 

5.2.  Clear  Water  Experiments 

We  first  apply  the  proposed  technique  to  the  clear  water  dataset  obtained  on  9/18/2013  [3].  Using  a  projection  path  beam 
concentration  optical  front  end  [3],  the  achievable  horizontal  image  resolution  was  512  pixels/line.  For  this  dataset,  the 
two  constants  in  Equation  (12)  were:  y=0.25  and  //= 0.25 .  The  comparison  of  the  performance  of  the  proposed  technique 
with  the  original  DCSa  based  approach  using  the  clear  water  dataset  are  shown  in  Figure  8.  The  proposed  technique 
achieves  two  improvements  over  the  original  DCSa  based  approach:  the  image  contrast  is  improved  due  to  the  reduction 
of  the  noise  in  the  black  regions,  and  a  better  image  resolution  is  apparent.  These  are  more  evident  by  comparing  the 
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high  resolution  regions  of  the  target,  marked  on  Figure  8a  and  8d.  Figure  8  illustrates  the  performance  at  several  different 
compression  ratios  (CRs).  The  image  contrast  and  resolution  improvement  can  be  observed  at  all  the  compression  ratios. 
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Group  =2 


(b)  128  meas/line,  CR=1:4 
Group  =4 
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Group  =2 


(e)  1 :4  compression,  Group 
=4 


(f)  1:8  compression,  Group 
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Figure  8  Comparison  of  the  performance  of  the  proposed  technique  with  the  original  DCSa  based  approach  using  the  clear  water 

dataset  at  different  compression  ratios 

The  AIQ  scores  in  Figure  9  reinforce  the  visual  observation  from  Figure  8.  Interestingly,  the  integrated  technique  produced  a  higher 
AIQ  score  at  CR=8:1  than  the  AIQ  score  at  CR=4:1  achieved  using  the  original  DCSa  technique. 


Figure  9  Comparison  of  the  AIQ  metric  of  the  proposed  technique  with  the  original  DCSa  based  approach  against  the  clear  water 

dataset  at  different  compression  ratios 

5.3.  Turbidity  cycle  tests 

We  also  apply  the  technique  to  the  turbidity  cycle  dataset  acquired  on  4/28/2014.  During  this  test,  a  different  illuminator 
design  that  focuses  the  beam  along  the  illumination  path  was  adopted  [3],  to  increase  the  system  power  output.  The  cost 
of  this  method  is  reduced  optical  front  end  robustness  against  external  vibrations  during  the  measurement  phase,  which 
yielded  the  more  noisy  regions  on  the  images  shown  in  the  following  figures.  The  achievable  horizontal  resolution  is  448 
pixels/line. 
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Figure  10  compares  the  results  of  the  two  techniques  at  C=0.6  (3  attenuation  lengths)  at  three  different  compression 
ratios.  The  two  constants  in  Equation  (12)  were:  y= 0.6  and  fi= 0.5. 
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Figure  10.  Comparison  of  the  performance  of  the  proposed  technique  with  the  original  DCSa  based  approach  using  the  test  dataset 

with  C=0.6,  3  attenuation  lengths,  at  different  compression  ratios 


The  proposed  technique  again  produced  both  contrast  and  resolution  improvements  at  different  compression  ratios.  The 
resolution  improvements  were  more  evident  in  the  marked  high-resolution  regions  in  Figure  10.  Such  improvement  is 
again  reinforced  from  the  corresponding  AIQ  scores  in  Figure  1 1  below. 


Figure  1 1  Comparison  of  the  AIQ  metric  of  the  proposed  technique  with  the  original  DCSa  based  approach  against  the  dataset  with 

C=0.6,  3  attenuation  lengths,  at  different  compression  ratios 


Figure  12  compares  the  results  of  the  two  techniques  at  C=0.8  (4  attenuation  lengths).  The  two  constants  in  Equation 
(12)  are:  y= 0.9  and  fi= 0.5.  As  with  the  case  of  C=0.6,  the  image  contrast  and  resolution  improvements  can  be  observed 
at  all  the  compression  ratios. 
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Group  =4 
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Group  =8 


Figure  12.  Comparison  of  the  performance  of  the  proposed  technique  with  the  original  DCSa  based  approach  at  different  compression 

ratios,  (C=0.8,  4  attenuation  lengths) 


The  AIQ  scores  in  Figure  13  below  confirm  the  visual  observations  based  on  Figure  12. 


Figure  13  Comparison  of  the  AIQ  metric  of  the  proposed  technique  with  the  original  DCSa  based  approach  against  the  dataset  with 

C=0.8,  4  attenuation  lengths,  at  different  compression  ratios 


6.  CONCLUSIONS 

This  paper  describes  a  new  reconstruction  technique  for  a  previously  proposed  compressive  line  sensing  underwater 
imaging  system  [1].  One  deficiency  in  the  existing  CLS  reconstruction  framework  is  that  the  interferences  from  the 
volume  backscatter  and  the  forward  scattering  from  the  adjacent  lines  are  not  adequately  mitigated.  The  core  idea  in  the 
proposed  technique  is  to  incorporate  the  Bayesian  inference  based  dynamic  compressive  sensing  approach  into  the 
original  “memoryless”  distributed  compressive  sensing  framework.  The  non-causal  reconstruction  framework  resulted 
from  such  integration  enables  us  to  better  mitigate  the  aforementioned  issues,  resulted  in  the  improved  image  contrast 
and  resolution.  The  effectiveness  of  the  technique  was  validated  using  data  acquired  in  the  experiments  conducted  in  the 
FAU  Harbor  Branch  optical  test  tank. 
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One  feature  of  the  CLS  imager  as  a  “coding  machine”  is  that  there  is  a  certain  level  of  disjoint  between  the  hardware- 
dependent  measurement  acquisition  (i.e.,  encoding)  process  and  the  image  reconstruction  (i.e.,  decoding)  process.  The 
work  described  here  reveals  a  benefit  of  such  separation.  With  the  same  data  acquisition  hardware,  it  is  possible  to 
enhance  the  system  performance  through  algorithm  improvement.  Such  adaptive  underlying  architecture  and  the 
algorithm-centric  implementation  can  potentially  reduce  innovation  and  enhancement  costs.  Although  a  more  robust 
prototype  system  should  certainly  improve  the  performance  as  well.  The  deficiencies  of  the  prototype  system  provided 
an  opportunity  to  demonstrate  the  ability  of  the  CLS  imager  to  “self-repair”  a  hardware  limitation  through  cost  function 
adaptation. 

The  proposed  technique  assumes  a  static  environment  (i.e.,  time  invariant  point  spread  function).  It  therefore  will  be 
interesting  to  study  how  to  further  improve  the  proposed  technique  so  that  it  can  be  more  effective  in  the  more 
challenging  cases  where  the  point  spread  function  is  time  variant  -  such  as  when  imaging  through  turbulence. 
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